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The finite-temperature optical conductivity in the planar t — J model is 

analysed using recently introduced numerical method based on the Lanczos diag- 
onalization of small systems (up to 20 sites), as well as by analytical approaches, 
including the method of frequency moments and the retraceable-path approxima- 
tion. Results for a dynamical mobility of a single hole at elevated temperatures 
T > t reveal a Gaussian-like i^(oj) spectra, however with a nonanalytical behavior at 
low uj. In the single hole response a difference between the ferromagnetic (/ = 0) 
and the antiferromagnetic (/ > 0) polaron shows up at T < /. At larger dopings 
numerical results in studied systems are consistent with the thermodynamical behav- 
ior for T > T* > 0.1 t. spectra show a non-Drude falloff at large frequencies. 
In particular for 'optimum' doping ~ 0.2 we obtain in the low-cj,r regime the 
relaxation rate ~ 0.6(a; -|- ^T) with ~ 3, being consistent with the marginal 
Fermi liquid concept and experiments. Within the same regime we reproduce the 
nearly linear variation of dc resistivity p with T. This behavior is weakly dependent 
on /, provided that J < t. 
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I. INTRODUCTION 



Unusual normal state properties of superconducting cuprates still lack a proper theo- 
retical explanation. There exists however a widespread consensus that most anomalies are 
closely related to strong electron correlations in these materials. The charge response, as 
manifested by the optical conductivity cr[ijj) and by the dc resistivity /?, has been extensively 
studied experimentally in cuprates [1,2], as well as theoretically (although with less conclu- 
sive results) within some models of strongly correlated systems. This paper is devoted to 
the study of the finite-temperature conductivity cr[ijj) in the prototype planar t — J model 
and represents an extension of the recently published analysis of the same problem [3]. 

The t — J model is one of the simplest models for strongly correlated electrons [4], 



where c-^{cis) are projected fermionic operators, taking into account that the double occu- 
pancy of sites is not allowed, rii and Si are corresponding local number and spin operators. 

Restricting our discussion to the two-dimensional and higher-dimensional systems, the 
most conclusive theoretical results have been so far obtained for the particular problem of 
a single mobile hole introduced into the reference insulator. The insulator is represented by 
the t — J model at half filling n = 1, i.e. with all sites singly occupied. In their seminal 
work Brinkman and Rice [5] solved the problem for the extreme correlation J = within 
the retraceable-path approximation (RPA). They pointed out on an essentially incoherent 
charge motion and evaluated dc mobility fio{T) of the hole. The latter exhibits a linear 
variation //q oc T, at least for high temperatures T > t. The same method has been applied, 
in connection with experimental results on cuprates, to the optical conductivity cr[ijj) [6]. 
Here, the incoherent motion shows up in the slow non-Drude falloff cr[ijj) oc l/cu for larger 
id. The RPA has been recently justified and applied more rigorously for infinite-D lattices 
[7]. An alternative approach is the evaluation of frequency moments of starting at 

r ^ t, first applied to the J = problem [U = oo Hubbard model) by Ohata and Kubo [8]. 
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On the other hand, the single-hole ground-state (T = 0) optical response has been in recent 
years well established by numerical studies of small systems via the Lanczos diagonalization 
method [9-11]. Nevertheless there are important unsolved questions, even for the single 
hole problem, e.g.: How accurate is the RPA for cr[ijj) on the planar lattice? Which are new 
qualitative spectral features at T < t, both for the J = and the J > case? 

More challenging open questions appear at finite hole doping. An evident one is related 
to the origin of the nearly linear variation p oc T in optimally doped cuprates [2]. Whereas 
at a low doping level one could possibly treat transport within the (semiconductor) model of 
independent quasiparticles, i.e. spin polarons, this concept could fail for moderate dopings, 
e.g. for Ufi > 0.1, due to the overlap between extended deformations around holes. An 
attractive proposal in the moderate doping regime is the one of spinous and holons as basic 
low-energy excitations [12], serving as the starting point for the gauge-theory approaches 
[13]. On the phenomenological level the marginal Fermi liquid (MFL) concept [14] has been 
very successful in linking many anomalous normal state properties of cuprates, in particular 
also cr[ijj) and p{T), by introducing in the Drude form u- and T-dependent relaxation rate 
T~^. Numerical studies of dynamical quantites within the ground state have been so far less 
conclusive in this regime [11]. 

The present authors introduced a numerical method [15], allowing for the calculation of 
dynamical response at T > 0. It enables more reliable studies of the normal state (higher-T) 
properties within correlated models, in particular due to a more controlled way of monitoring 
finite size effects. The method has been applied to the finite-T optical conductivity (T{io) 
[3], and more recently also to the spin response [16] in the t — J model. Some important 
conclusions, in particular those concerning the comparison with experimental results in 
cuprates, have been given in the short communication [3]. We present here a more complete 
numerical analysis of the optical and the dc conductivity. 

The organization of the paper is as follows: In Sec. II we calculate analytically several 
frequency moments (up to the fourth order) of (T{io) at T ^ t within the t — J model on 
a square lattice. For comparison we also analyse the single hole dynamical mobility within 
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the RPA. Sec. Ill is devoted to the presentation of the numerical method and related tests. 
In Sec. IV we first analyse the dynamical and the dc mobility of a single hole, both for J = 
and J > 0. Next the results for cr[ijj) and p{T) at finite hole doping are considered. In 
particular we extract the phenomenological relaxation time t{lo,T), which we compare to 
the MFL theory. 



II. ANALYTICAL APPROACHES 



A. Moment method 



Within the linear response theory the optical conductivity (T{io) is given by 

a{u) = CH, C{co) = elRe dr e'^^ {j {r) j) , (2) 

to Jo 

where /3 = l/ksT and j is the particle current density along the chosen axis, e.g. the x axis. 
In the t — J model (1) with N sites, j can be expressed as (we use in the following units 
with h = ks = 1 &s well as the lattice spacing ao = 1), 

tt 



J 



The method of frequency moments has been one of the standard ones to analyse the 
dynamical response [17]. It has been used for the conduction in a correlated system first 
by Ohata and Kubo [8], where cr[ijj) in a simple cubic lattice has been analysed. Analogous 
approach, combined with the high-T expansion, has been recently used to evaluate other 
dynamical quantities in the t — J model [18]. We are using in the following the method for 
the t — J model on a square lattice and restrict our analysis to high temperatures T ^ t. 

The goal is to calculate frequency moments of 



= dto to "(7(0;) = 7r/3eoa2n , (4) 

00 

taking into account that odd moments vanish. For T ^ t moments are related to the 
short-time expansion of the current-current correlation function, i.e. 



{Jir)j) = E Vl"^'^^"^ ' 

«2. = Tr(j('^)j('^))/Tr 1, j^^^ = H^j, (5) 

and HA = [H,A]. 

It is instructive to obtain analytical expressions for few lowest moments since they 
provide an insight into the processes leading to the current relaxation. The evaluation of the 
moments is straightforward, but very tedious for larger n. We follow a procedure analogous 
to the one used in [8], summing (tracing) in Eq. (5) over all initial configurations and over 
processes which restore these configurations. 

The zeroth order moment is independent of the lattice dimensionality (in hypercubic 
lattices) [8]. It requires a simple hopping process, an occupied site and an unoccupied site 
being neighbors, 

ao = Tr(jj)/Tr 1 = 2fnenh, (6) 

where we introduce Ue = NpjN and nh = l — np = N^/N as the electron and the hole 
density (per site), respectively. 

The second order moment involves contributions from the hopping and the spin part in 
Eq. (1). Let us first consider the hopping contribution, where we can represent both Hf and 
j as a sum over hoppings between directed pairs of sites, 

Ht = J2 H^J , J =iJ2 e^JH^J, (7) 

and 

jI'^ = [H,,j] = H,j-jH,. (8) 

It is important to realize that due to Eq. (7) vanishes identically in the ID chain of corre- 
lated electrons [5], as well as for hoppings along one direction only, being a manifestation of 
the current conservation in ID. On a square lattice nonvanishing contributions to come 
from the processes on a plaquette, as represented in Eigs. 1(a), (b), under the condition that 
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fermionic spins on the plaquette are not equal. Summing over different incoherent contribu- 
tions to Eq. (5) within a given configuration and different realizations of the plaquette, we 
obtain 

a2t = Gt^nlrih + St^nlnl, (9) 

where the first and the second term come from the configurations as shown on Fig. f(a) 
and Fig. f(b), respectively. We treat the exchange term Hj in an analogous way and one 
contributing process to j^p is shown in Fig. f (c). The result is 

«2J = 4tV2 J2 iiS^ ■ Sif) = -/J'nlnu . (10) 

/ n.n. i 

The fourth order moment requires quite tedious calculations, hence we restrict ourselves 
to particular cases: either the vicinity of half-filling n/^ ^ 1 or to nearly empty band ne ^ 1, 
both in the strong correlation limit J = 0. We have to consider the action of the operator 

= [H,,j[^\ = H,H,j - 2H,jH, + jH,H, . (11) 

Again, hopping contributions along only one chosen axis vanish due to = 0, leading 
to a4 = in ID. The main task is to group all processes which could lead to the same 
intermediate state through the operation of j^^^ and thus contribute coherently to a^. Two 
examples for Uh <^ 1 are in Figs. 2(a), (b). The coherence factor depends on the particu- 
lar spin configuration. F.g., whereas in Fig. 2(a) three different final states are obtained 
(resulting in the coherence factor C = 6), in Fig. 2(b) only two final states are different, 
leading to a partial cancellation and C = 2. The enumeration of all nonequivalent diagrams 
(analogous to those in Figs. 2(a), (b)) leads to 

04 = QQt^rih , J = 0, nh<.l. (12) 

An analogous consideration of processes involving configurations with two fermions (with 
opposite spins analogous to Fig. 1(b)) on a rectangular plaquette yields 

04 = iSt^nl , J = 0, ne < 1. (13) 
6 



What can one learn Irom moments, as presented above lor T ^ t 1 At low hole doping 
Ufi <^ 1 one can interpret conductivity in terms ol the dynamical mobility fi{io) or diffusion 
response -D(cu) ol independent holes, related by 

cr{ijj) = eonhjJ,{io), jJ,{io) = /3eo-D(cu), (14) 

where the latter relation is meaninlul only lor T ^ t. The diffusive motion ol holes in 2D (in 
contrast to the nondiffusive one in ID) emerges due to the noncommuting processes around 
individual plaquettes, as presented in Figs. 1,2. The lowest frequency moments ol -D(cu) lor 
Hh < 1 follow Irom Eqs. (6,9,10,12), 

ao 16 

(u*) = ^ = 33t^ J = . (15) 

From Fq. (15) we note that the exchange term contributes less to the hole scattering (pro- 
vided that J < t) than the strong correlation requirement ol no double occupancy. It is also 
interesting that lowest moments predict a nearly Gaussian form ol -D(cu), since (u'^) ~ 3(cu^)^. 
One can therefore extract an approximative dc diffusion constant Dq as 

D. ~ ^£L^ = f = 1.447f , (16) 

to be compared later with other methods. 

At low electron density ne ^ 1 the current scattering is essentially different. The Ire- 
quency moments 

{lo^) = 4:t\,, {uj^) = 24t^ne (17) 

are non-Gaussian (cu'*) ^ 3(c<;^)^. It is easy to recognize that processes included in Fqs. (9,13) 
contain only two-body collisions ol electrons with opposite spins. For configurations with a 
pair ol electrons with parallel spins the current can decay only due to the processes involving 
other electrons, hence terms with a4 oc are crucial etc. Therefore, one can expect quite 
different spectra cr[ijj) with possibly more singular dependence at io <^ t. 
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B. Retraceable path approximation 



The RPA has been introduced [5] for the analysis of a single-hole motion in the strongly 
correlated system with J = 0. Assuming that the motion of a hole in the spin background 
is entirely incoherent, i.e. without closed loops, only retraced paths contribute to the local 
propagator G'(cu), approximated as 

G{u) = , (18) 

[z — 2)ijj + zJio'^ — ujq 

where z is the coordination number of a lattice (on a square lattice z = A) and loq = 



z — 1 t. The RPA is exact in ID, on Bethe lattices, as well as on infinite-D lattices 
[7]. Within the RPA it is possible [5,6] also to evaluate the optical conductivity c(cu). The 
mobility, defined by Eq. (14), can be expressed as 

fi{io) = 2f- / duj'e-^'^yiuj^uj + uj')/ duj'e-^''G"iuj'), (19) 

(jJ J —iUQ J —iUQ 

where 

J^{(jJl,(jJ2) = ^ SiS2F{(jJi + iSi6,(jJ2 + is2S), 

Sl ,S2=±1 

F{co,,co2) = GMG{co2)[l + — "'^"^"^.J , (20) 

aia2 aia2(c(iC(2 ~ tj 

«1,2 = ^[^1,2 + \/^l,2-^o] ■ 

For T ^ t the diffusivity spectra -D(cu) = Tfi[ijj)/eo , resulting from Eqs. (18-20), are 
presented together with numerical result in Fig. 6. In spite of the finite cutoff at cOmax = 
2ijJo = \fY2 t, spectra have a close resemblance to the Gaussian function. This is clear from 
the lowest moments within the RPA, 

iu?) = u\ (u*) = m\ (21) 

which satisfy the Gaussian relation, but are as well close to the exact ones, Fq. (15). The 
agreement is reduced if one uses in Fq. (20) the simplified form without the correction factor, 
i.e. F[ijJi^ijj2) = G{iLii)G{iLi2). The latter approximation yields (cu^) = 8t^, (cu'*) = 152^^. 
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We also note that D[ijj = 0) = 1.367 t within the RPA, in comparison to 1.447 t, obtained 
by the moment method, Eq. (16). 

III. NUMERICAL METHOD AND TESTS 

An efficient numerical method which allows the study of T > static and dynamical 
quantities in small correlated systems, has been recently introduced by the present authors 
[15], and applied to the dynamical conductivity [3] and to the spin susceptibility [16] in the 
t — J model. The method, here used for the evaluation of the current-current correlation 
function C(cu), is based on two essential ingredients: 

(a) the calculation of approximate eigenfunctions l?/;^) and IV'fc), A; = 1, Af, with correspond- 
ing energies €nk and e„/;, generated by the Lanczos diagonalization procedure from the initial 
functions |n) and respectively, and 

(b) on a reduction of the full trace to a partial random sampling over basis functions 

No M 

C{u) = ^-^7re^^^e-'^-'"(n|C)(Cb#n(V5fc1j»% - e.™ + ink). 

n m,k 

No M 

z = j:j:--''-\Hrj\'. (22) 

n m 

It is important to note that both A^o and M are much smaller than the full-basis dimen- 
sion Nfot. hence the requirements for the T > calculation are comparable to the T = 
Lanczos method for response functions [9-11]. An important feature of the method is that 
it reproduces correct frequency moments of cr[ijj) (for a finite system) at higher T as well as 
a meaningful T = result. 

Some tests of the method, in particular in the connection with the problem of dynamical 
conductivity, have been presented already in Refs. [15,3]. When discussing the feasibility 
of the method two different aspects have to be considered: (a) the role of approximations 
within the chosen finite-size systems, and (b) the finite size effects due to the smallness of 
the systems studied. 
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(a) Approximations introduced by using Af, A^o ^ ^tot seem to be quite well understood 
and controlled. At full sampling A^o = ^tot the number of Lanczos steps M determines 
the number of correct frequency moments (at least for higher T) within the chosen system, 
or equivalently the frequency resolution in spectra to Aco ~ cOmax/M^ ^max representing 
typically the energy span of the model. Since the information content in these moments is 
limited, in particular due to finite size effects, there is no point in using very large Af, i.e. 
we restrict our calculation in most cases to Af < 240. The effects of using reduced sampling 
A^o ^ A^tot are most pronounced at low T, where only symmetry sectors with the lowest 
energies contribute. 

We test the performance of the method on the example of a system of 10 sites and 
Nh = 3 holes. While the full basis dimension is A^tot = 15360 in this case, we choose Af = 40 
and A^o = 180 in our method, i.e. M^Nq <^ Ntot- We compare results to the exact ones, 
calculated via the full diagonalization, which requires the diagonalization of the 420 X 420 
matrix in the largest Sz = ^f^^q = sector. In Fig. 3 we present cr[ijj) spectra at different 
r/t, obtained by both methods (plotted with an additional broadening of rj = 0.05 t). It 
should be noted that the smaller sampling A^o = 180 represents within the ground state 
sector only A^i = 45 states. In spite of a small effective number of samples even the T = 
result is (surprisingly) robust, indicating that the response of different samples is quite 
similar. It is also plausible that the effects of finite sampling are even less pronounced at 
higher T, as evident from Fig. 3. 

(b) Turning to finite size effects an important conclusion can be drawn from Fig. 3. 
Fven at r > t, spectra for the 10-site system exhibit nonuniversal features, which could be 
attributed to finite size effects. The latter entirely disappear in larger systems, e.g. for A^ = 
16 (at same T) as presented in Sec. IV. The difference is clearly due to the dramatic increase 
in the number of states Ntot with A^, whereby the number of peaks in the spectra scales even 
as N^^^ (within our method as A^o^^)- This leads to a large reduction of (finite-size induced) 
fiuctuations in spectra at particular cu, requiring very small additional broadenings rj <^ J 
to average them. Thus one of the crucial advantages of the new method is in its ability to 

10 



study systems which yield very smooth, macroscopic-like dynamical spectra at T > 0. 

The density of (many-body) states becomes very sparse on approaching the ground state. 
Consequently, the thermodynamic averaging becomes ineffective at T < T* oc AEq^ where 
AEq is an average level spacing in the low-energy sector. Generally, the crossover at T* can 
be monitored by the emergence of finite size gaps in spectra etc. For cr[ijj) an additional test 
is the charge stiffness [19], 

1 1 

= -^j{HMn,xx) o / Cr{to)dLO, (23) 

ZiV TTCq Jo+ 

which measures the coherent charge propagation around the system. For a macroscopic 
system at T > (with the current being a nonconserved quantity) one expects generally 
-De = 0, although this does not appear to be the case for (nontrivial) integrable quantum 
many-body systems [20]. 

In Fig. 4 we present results for Dc{T) for a single hole in the t— (J = 0) model. This 
example is particularly simple, since the ground state is known to be ferromagnetic (FM) 
with Ekin = — 4t as the consequence of the Nagaoka theorem [21]. Hence at T = only 
the coherent contribution at = remains in c(cu), and Dc{T = 0) = t/N. Fig. 4 shows 
that Dc{T) interpolates quite smoothly between Dc = at T = oo and Dc{T = 0). The 
crossover however becomes increasingly sharp at larger systems > 18, where it is impor- 
tant to notice the qualitative improvement when doubling the size from = 10 (studied 
mainly so far by the method of full diagonalization [22]) to = 20. Fven for the largest 
system = 20 we find that for T < 0.2t cr[ijj) becomes dominated by the undamped co- 
herent contribution. Therefore, in this regime low-cu spectra do not represent properly the 
thermodynamic behavior. 

For J = 0.3 t the results for Dc{T) are less regular. We show in Fig. 5 results for various 
Nfi > 1 on a system of A^ = 16 sites. When judging the extent of deviations of Dc from zero 
at r > 0, it is useful to compare values with the maximum possible ones, i.e. with the value 
of the sum rule at T = 0, D^ax = \Ekin{T = 0)|/4A^, as follows from Fq. (23). We notice 
from Fig. 5 that Dc typically shows a rather abrupt transition from Dc ~ (small deviations 
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could be attributed to finite sampling and to finite width of cu-channels) to Dc{T) ^ at 
the crossover temperature T*, depending on the hole doping in the particular system. 
For T < T* the variation Dc{T) can become quite unphysical, i.e. we get in some cases even 
D < 0. These phenomena are infiuenced by the particular use of boundary conditions and 
more sensible results can be obtained by the introduction of additional fiuxes (or twisted 
boundary conditions) [19]. Nevertheless we are here interested only in the regime T > T*. 
From Fig. 5 it follows quite clearly that T* is minimum, i.e. T* ~ 0.1 t, for the intermediate 
(optimum) doping regime 0.1 < < 0.3, and is larger for the underdoped as well as for the 
overdoped systems. 

IV. NUMERICAL RESULTS 
A. Single hole mobility 

The properties of a single hole are the basis for the understanding of the conductivity in 
the low-doping regime. Let us first consider results for the prototype, but notrivial, problem 
of the dynamical diffusion -D(cu), Fq. (14), at J = and T ^ t. The spectrum, as presented 
in Fig. 6, has been obtained for = 20 and M = 200 to achieve high frequency resolution. 
For comparison we show also the RPA result, calculated from Fqs. (19,20). 

It should be noted that for = 20 we reproduce already the exact frequency moments 
(for an infinite system) up to the fourth order, Fq. (15), i.e. (cu^) = and (cu'*) = 33^^. 
The overall agreement of the RPA result with the numerical one is reasonable, taking into 
account the approximations involved in the RPA. It is however quite evident from Fig. 6 
that in contrast to the RPA the actual -D(cu) is nonanalytical at io = 0. We reproduce the 
same phenomenon in A^ = 18 system, but not within the same systems with J > 0. We do 
not have a proper explanation for this effect. A possible interpretation can be given in terms 
of finite ferromagnetic-like clusters. Probability for larger FM regions is decaying fast, on 
the other hand holes exhibit a nondiffusive motion in these regions, so that the interplay of 
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both effects coufd resuft in a nonanaiyticai D[ijj ~ 0). Such behaviour is not uncommon in 
the transport phenomena, Icnown e.g. in connection with iong-time taiis [23]. 

The spectra for J = 0, but finite T < oo, are presented in Fig. 7(a). Piotted are 
dynamicai mobiiities Eq. (14), again caicuiated for N = 20. Comparing Fig. 7(a) 

with Fig. 6 we conclude that for T > t the spectra essentially retain the high-T form, i.e. 
fi{io) ~ D{lo,T = oo) jT. Spectra change qualitatively for T < where the quasiparticle 
central peak starts to emerge at low cu, and the incoherent broad background (peaked at 
~ t) vanishes on approaching T — )• 0. In this case the relevant quasiparticle at T ^ t 
represents the FM polaron. However, our systems are too small to be able to obtain some 
more quantitative conclusions on the FM polaron relaxation. 

J > does not infiuence spectra essentially at T > t, since according to Fq. (15) 
J < t only slightly increases the width (cu^). As seen from Fig. 7(b), there are qualitative 
differences appearing for T < J. Namely, for T — )• we are dealing with a nontrivial 
ground state of the antiferromagnetic (AFM) polaron, which has an enhanced mass but a 
nonvanishing incoherent spectral component at io > J. 

Finally we extract from dynamical quantitites the dc values //q = /w(^ — > 0), or equiva- 
lently an effective diffusion constant Dq = T/io/eo (which is not equal to a proper diffusion 
constant for T < t). It is plausible that the dc values can be obtained in a unique way 
only for T > T* where Dc <^ D^ax- In Fig- 8 we present numerical results for DoiT), both 
for J = and J = 0.3 t, whereas the RPA (J = 0) result is shown for comparison. It is 
plausible that at T > t finite J > slightly reduces Dq, as follows also from Fqs. (15,16) 
applied to J > 0. On the other hand, both numerical and the RPA result show a decreasing 
tendency of Do[T) with decreasing T < t. In the resistivity p{T) oc T / Do[T) this would 
show up as the deviation from the linear T-variation, i.e. p{T) would remain rather constant 
(in particular for J = 0.3 t) in the intermediate regime T* < T < t. At low temperatures 
T < t the relative reduction of Dq due to J > is much larger, indicating that the scattering 
of the hole on AFM excitations (magnons) is more efficient. 
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B. Conductivity at finite hole dopings 



Results lor a single hole = 1 can be interpreted also as low-doping results for the 
(sheet) conductivity cr[ijj) in finite systems. In this sense we presented in [3] spectra for 
Ufi = 1/18, which are essentially equivalent to the results for = 1/20 shown above. By 
increasing the doping, cr[ijj) could reveal new features, even at T > T* . It is convenient to 
discuss in this case the dimensionless quantity ct/Jq, where po = fij = 4.1A;0 is the universal 
2D sheet resistance. 

We first show in Figs. 9 optical conductivities at fixed J = 0.3 t, but at various dopings 
Tih = 2/16,3/16,8/16. Related spectra for rih = 4/16 have been displayed already in [3], 
being very similar to the rih = 3/16 case. One feature common to all rih considered (including 
Nh = 1) is the non-Drude type behavior for u > t. This confirms the belief that the 
incoherent motion of holes (dominating the high-cu response), as e.g. described within the 
RPA [6], can remain a valid concept even at large hole doping rih ~ 0.5. 

Some qualitative differences appear at low u. Results for = 1 indicate that for 
T < T* ^ 0.7 J a quasiparticle peak and a related quasigap starts to emerge at low cu, as 
seen in Fig. 7(b). Some remains of this phenomenon could persist even at = 2/16 in 
Fig. 9(a). On the other hand, no such effect is observed in the regime of 'optimum' doping 
with Ufi = 3/16,4/16, where no additional low-cu scale can be resolved down to the lowest 
T = T* ^ 0.1 t. On increasing further, T* increases substantially, indicating that the 
scattering becomes again less effective and the mean free path (at fixed T) longer. Hence, 
we cannot efficiently investigate low-electron densities Ue < 0.5, where one could presumably 
expect the Landau- Fermi-liquid behavior at low frequencies. 

We calculate also the corresponding spectra for J = 0. Two of them are presented in 
Figs. 10. It is plausible that for T > J spectra are not significantly infiuenced by J. Also 
at high dopings > 0.4, J becomes irrelevant at arbitrary T (provided that J < t). F.g., 
for Tih = 8/16 we do not observe any significant difference between the J = 0.3 t and the 
J = spectra. There appear however some systematic changes at T < J within the low 
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doping regime < 0.2. Central peaks at the origin in Figs. 10(a), (b) (at T < 0.5t) seem 
to be sharper and more distinguishable than the ones in Figs. 9(a), (b). One could possibly 
attribute this difference to the known tendency towards the FM (at least partially) polarized 
ground state at J = 0, persisting even at > 0. 

With respect to experiments on cuprates [1], the most challenging is the case of 'optimum' 
hole doping ~ 0.2 within the AFM correlated spin background generated by J > 0. In 
this regime optical measurements reveal a very universal behavior of c(cu). In the previous 
paper [3] we showed that our model results for = 4/16 reproduce (assuming J/t = 
0.3, t = 0.4 eV) quantitatively sheet conductivities in BISCCO, YBCO [24] and LSCO 
[x = 0.2) [25]. It is very instructive to analyse the complex cr[ijj) in terms of the effective 
frequency dependent relaxation time t(c<;) and the effective mass m*(c<;), introduced via the 
memory function M [26] , 

^^""^ = u + M{u) ' ^ = -{HHn,..)/N, (24) 

and 



M"{ 



-{uj) l + M'{uj)/uj' 



(25) 



nit S ' 

where nit = l/2taQ is the bare band mass. Using Fq. (25) one can formally rewrite (T{io) in 
the familiar Drude form 

Using Fq. (25) we can evaluate from known cr[ijj) both r and ni*. Results for = 
3/16, J/t = 0.3 are presented in Fig. 11. From Fig. 11(a) it follows that in the whole 
low-cu, T regime we can remarkably well describe the behavior of r with a simple law 

= 27rA(cu + Cr), A -0.1, C~3. (27) 
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This is just one of the forms proposed within the MFL theory [14]. It should be also noted 
that obtained A is very close to the experimental one A ~ 0.11 [27,1]. 

Also m* /nit in Fig. 11(b) qualitatively agrees with experimental findings [24,27,1]. For 
quantitative comparison we should note that by lowering T to experimentally investigated 
ones we expect the increase of m* at low io. On the other hand, assuming realistic values 
ao = 0.38 nm and t = 0.4 eV [4] we should take into account that nit ~ 0.6 niQ. 

Finally, we display in Fig. 12 results for dc resistivities p{T), as extracted from (t{lo — )• 0) 
(mainly presented and discussed already in [3]). We note that for T > t and < 0.5 
slopes are approximately dp/dT ~ CPo^b I n-ht with C, ~ 0.4, in comparison to the RPA 
value C ~ 0.72, as deduced from Fig. 6. It is remarkable that p{T) results are only weakly 
infiuenced by J down to T ~ J, whereby it is plausible that the effect of J nearly vanishes 
on approaching the electron-like regime < 0.5. 

V. CONCLUSIONS 

Main conclusions of numerical and analytical investigations, as presented above and in 
the preceding paper [3], can be summarized as follows: 

(a) The finite-T method yields within small systems meaningful (i.e. without pronounced 
finite size effects) spectra (T{io) for T > T*{nh). For the conductivity problem it is plausible 
to identify T* with the temperature, where the transport mean free path /q reaches (exceeds) 
the size of the system, i.e. Iq{T*) ~ \/iV- This fact represents an independent information 
on transport. Within the studied system of 16 sites it is evident that T* depends quite 
significantly on doping. It is lowest for the 'optimum' doping rih ~ 0.2 where T* ~ 0.3 J. 
In the low doping regime rih < 0.1 we find T* ~ 0.7 J, while in the overdoped cases with 
Tih ~ 0.5 (where J is irrelevant) T* ~ 0.5 t. 

(b) The single-hole dynamical diffusivity D{lo) at large T ^ evaluated from the lowest 
frequency moments or within the RPA, appears to be close to the Gaussian form. Numerical 
results confirm general features of analytical approximations, but for J = also indicate on 
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a fundamental nonanalytical behavior of D[ijj — )• 0). 

(c) The single hole dynamics becomes quite sensitive on J at low T < i.e. the hole 
can behave either as FM polaron for J ~ 0, or as an AFM polaron for J > 0. 

(d) At finite doping, beyond the low-doping regime rih > 0.15, both cr[ijj) and p{T) do 
not reveal an essential difference between J = and J = 0.3 t, at least not for T > T*. 

(e) At 'optimum' doping ~ 0.2 cr[ijj) spectra show a very universal form. For this case 
we extract the phenomenological relaxation time r and the effective mass m*. The form of 
t(c<;, r) is consistent with the MFL proposal and moreover agrees with experimental values. 

(f) Results for cr[ijj) and p{T) agree well with experiments without any adjustable pa- 
rameters [3] (both J and t are well fixed either by experiments or by other considerations). 
This confirms the belief that the strong correlation effects, as embodied e.g. within the t — J 
model, are crucial and possibly sufficient for the understanding of normal state properties 
of the cuprates and related materials. 
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Figure captions 

FIG. 1. Some processes with nonvanishing contributions to j^^K Initial configurations 
are represented as: empty sites by open dots, fermions with up spins by full dots, and 
fermions with down spins by crossed dots. Straight lines show hopping and wavy lines spin 
flip processes. Numbering represents the order of processes. 

FIG. 2. Some processes with nonvanishing contributions to j^'^K 

FIG. 3. cr[ijj) spectra for = 3 holes on = 10 sites obtained via the full exact 
diagonalization (full line), and by using our method with the sampling A^o = 180 and 
Af = 40 Lanczos steps (dashed line). 

FIG. 4. Charge stiffness Dc{T) for a single hole in the t — (J = 0) model for planar 
systems of various sizes. 

FIG. 5. Normalized charge stiffness D^ax for a various number of holes in the 
model with 16 sites and J jt = 0.3. 

FIG. 6. Dynamical diffusion D{lo) for a single hole at T = oo and J = 0, as obtained by 
the numerical calculation on a 20-sites system (full line), and the RPA theory (dashed line). 

FIG. 7. Dynamical mobility /u(cu) for a single hole at various T in the planar t — J model 
with: (a) J = 0, and (b) J = 0.3 t. 

FIG. 8. Fffective diffusion constant Do[T) for a single hole. Shown are numerical results 
obtained on a system of = 20 sites for J = (circles) and J = 0.3 t (squares). Dashed 
line represents the RPA result. 

FIG. 9. cr[ijj) spectra for J = 0.3 t and various T at finite hole dopings: (a) = 2/16, 
(b) Ufi = 3/16, and (c) = 8/16. The smoothing width here and in Fig. 10 is ?y = 0.07 t. 

FIG. 10. cr[ijj) spectra for J = and various T at finite hole dopings: (a) = 2/16 and 
(b) Hh = 3/16. 

FIG. 11. (a) Fffective inverse relaxation time t(c<;), and (b) the mass enhancement 
m*{Lo) I nit, for various T in the system with rih = 3/16 and J = 0.3 t. 
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FIG. 12. Resistivities p{T) at different dopings for J = 0.3 t (full lines) and J = 
(dashed lines). 
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